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We present a detailed study of the finite temperature dynamical properties of the quantum Potts 
model in one dimension. Quasiparticle excitations in this model have internal quantum numbers, 
and their scattering matrix deep in the gapped phases is shown to take a simple exchange form 
in the perturbative regimes. The finite temperature correlation functions in the quantum critical 
regime are determined using conformal invariance, while far from the quantum critical point we 
compute the decay functions analytically within a semiclassical approach of Sachdev and Damle [K. 
Damle and S. Sachdev, Phys. Rev. B 57, 8307 (1998)]. As a consequence, decay functions exhibit 
a diffusive character. We also provide robust arguments that our semiclassical analysis carries over 
to very low temperatures even in the vicinity of the quantum phase transition. Our results are also 
relevant for quantum rotor models, antiferromagnetic chains, and some spin ladder systems. 



PACS numbers: 05.30.-d,05.50.+q,73.43.Nq 
I. INTRODUCTION 

Quantum critical systems have been in the focus of in- 
tense studies in the last two decadesi^ In these systems a 
second order phase transition driven by quantum fluctua- 
tions takes place at T = temperature, which, neverthe- 
less, determines the physical properties of the system over 
a finite parameter range, and leads to anomalous dynam- 
ical scaling. Much of our interest in these quantum criti- 
cal systems has been motivated by the anomalous scaling 
behavior that can be observed in various heavy fermion 
compounds but the anomalous behavior of the nor- 
mal states of cuprates has also been interpreted in terms 
of an underlying quantum phase transition, possibly hid- 
den by the d-wave superconducting state, similar to many 
heavy fermion compoundsi&S Much of the studies of 
quantum critical phase transitions focused on the sim- 
plest case of magnetic phase transitionsjiaiSiii*^ where 
typically quantum fluctuations compete with a tendency 
of magnetic ordering. While in many cases a quantum- 
classical mapping enables one to map out the phase di- 
agram and determine critical exponents based upon our 
knowledge and experience with classical criticality, com- 
puting real time response functions represents a major 
challenge for theorists. In this regard 1+1-dimensional 
systems are of crucial importance since powerful methods 
can be used there to analyze their dynamical properties 
both in the gapped phases and in the quantum 
critical regimei. These one-dimensional models, besides 
being relevant to some experimental systems^ serve also 
as test grounds for higher-dimensional systems^ 

Maybe the simplest and most thoroughly studied 1+1- 
dimensional model is the transverse field Ising model, 
where a one-dimensional Ising chain with ferromagnetic 
coupling J is coupled to a magnetic field h in the x 
direction. In this model a quantum phase transition 
takes place from a ferromagnetically ordered to a para- 
magnetic state at a critical value of the magnetic field, 



h c /J = 1, and the critical theory is simply that of 
the two-dimensional classical Ising modcliiii^ In this 
model, dynamical properties arc well-understood both in 
the critical regime and in the gapped phascsii^ 

In case of the transverse field Ising model, our deep un- 
derstanding of the quantum critical regime is based on 
the observation that a 1+1-dimcnsional critical system is 
conformally invariant, and one can therefore describe it 
by means of some conformal field theory^ More interest- 
ingly, however, we can turn this observation backwards 
by saying that any conformal field theory should cor- 
respond to some model that displays a quantum phase 
transition. 

In this way, we can, in principle classify, construct, and 
study new one-dimensional quantum mechanical systems 
which belong to different universality classes, and thus 
exhibit new and interesting critical behavior. 

Relying on this observation, in the present paper we 
shall study the simplest possible 1+1-dimensional gen- 
eralization of the transverse field Ising model which be- 
longs to a different universality class, the Q-state quan- 
tum Potts model. In this one-dimensional model each 
spin has Q different components corresponding to the 
basis states \n) with /i = 1, . . . , Q, and the Hamiltonian 
takes on the following simple form 



(i) 



where the operator P M projects on state |/Lt) within the 
Q-dimcnsional local Hilbcrt space, while P projects to 
the 'corner' state ^ p ItA/VQ- The ferromagnetic cou- 
pling j thus tries to polarize all Potts spins at each site 
to one of the Q possible orthogonal directions, while 
the magnetic field h = jg tries to project each spin to 
the (1,1,...) direction and thus to destroy long range 
order generated by the ferromagnetic coupling. In the 
particular case of Q = 2 Eq. (|TJ reduces to the Ising 
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Hamiltonian, with £ M =i,2 Iflf+i = («+i + l)/2 and 
Pi = + For later convenience, we shift the 

ground state energy and replace the projectors in Eq. 
by the traceless operators P M = P M 1 — -^ and P = P — -^ . 

The quantum Potts model defined by Eq. Q and 
the corresponding two-dimensional classical Potts model 
have been studied extensively before ^i^iSSiSiiSSiS How- 
ever, while we know a lot about the quantum Potts 
model's thermodynamic behavior, much less is known 
about its real time dynamics, especially at finite temper- 
atures. The purpose of the present paper is to study in 
detail dynamical correlations within the quantum Potts 
model at finite temperatures, both in the gapped and in 
the quantum critical regimes. As we shall see, many of 
our results also apply to other one-dimensional gapped 
systems and are therefore of direct physical relevance. 

The one-dimensional quantum Potts model, Eq. Q, is 
known to exhibit a phase transition as a function of the 
coupling g02£&: Below a critical value, g < g c the ground 
state of Eq. JTJ is a Q-fold degenerate ferromagnet corre- 
sponding to the Q different ferromagnetic alignments of 
the Potts spins. The elementary excitations are domain 
walls that move along the chain with a dispersion 



e(k) 



(2) 



where the two indices \i and \j! denote the orientations 
of the ferromagnetic order parameters on the two sides 
of the excitation, and k is its lattice momentum. These 
quasiparticles are gapped and their energy can be ap- 
proximated for small g-s as e(k) ~ j — g ^ cos kj . 

For g > g c , on the other hand, the ground state is a 
non-degenerate paramagnet, which corresponds to orient- 
ing all spins along the (1, .. . , 1) direction. Elementary 
excitations of this state consist of A = 1, . . . , Q — 1 pos- 
sible local 'spin flips', which propagate along the chain 
with a dispersion 



(3) 



These excitations are also gapped and for very large val- 
ues of g one finds e(k) ~ j 9(1 — ^- cos A;) from pertur- 
bation theory. 

Clearly, since the structure of the ground state and 
its elementary excitations is entirely different for j< 1 
and j » 1, a phase transition must occur at some crit- 
ical value g = g c . This phase transition turns out to be 
of second order for Q < 4, while it is of first order for 
Q > l 20 i 24 . This implies that the quasiparticle gap re- 
mains finite for Q > 4, while it approaches zero at g = g c 
for Q < 4. The case Q = 4 is special, and shall not be 
discussed here: Although the phase transition is contin- 
uous for Q = 4 too, there a marginal operator dominates 
the critical points The case Q = 3, on the other hand, 
is of special interest: For Q = 3 one finds a standard sec- 
ond order quantum phase transition, which belongs to a 
universality class different from the Q = 2 Ising casein. 
Furthermore, for Q > 3 the quasiparticles have internal 



quantum numbers making their dynamics very interest- 
ing, and in many ways similar to that of quantum rotor 
models and S = 1 spin antifcrromagnetic chain o 1 ! 26 . 

The phase diagram of the Q = 3 quantum Potts model 
is sketched in Fig. ^ 
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FIG. 1: (Color online) Sketch of the phase diagram of the one- 
dimensional Q = 3 state quantum Potts model. The quantum 
fluctuations become dominant at a temperature scale T ~ A. 

As we indicated in the figure, the quasiparticle gap A 
approaches zero as a power law at the quantum critical 
point g c , 



A Q=3 (g) ~C±\g-g, 



|5/6 



(4) 



where C± denote prefactors corresponding to the g > g c 
and g < g c regimes. The exponent in Eq. (J3J sim- 
ply follows from the quantum classical mapping between 
the two-dimensional classical Potts model and the one- 
dimensional quantum Potts model)2i2i and is simply the 
critical exponent v = 5/6 of the correlation length £ ~ 
1/A in the two-dimensional classical problemiiiiSiSSiS In 
the quantum critical regime, A < T < j, physical prop- 
erties of the model are governed by the zero temperature 
critical point g = g c . 

The basic difference between the above two quantum 
phases can be captured by the change in the dynamical 
structure factor {ui, q), defined as the Fourier trans- 
form of the spin-spin correlation function, 



^^ty/J,' if: "^i) 



Pt(t) P$ 



(0) 



(5) 



The structure factor S*^ (to,q) is directly measured by 
neutron scattering in magnetic systems, and it is related 
to the dynamical susceptibility, 



Xwfrxi) = H) (tfTM ^ (0)]) 6{t) , (6) 
through the relation 



Sw>(u,q) = -2(n(w) + 1) Im{x M)M / (w, q)} 



(7) 



with n(u>) the Bose function. 

In the ferromagnetic phase at T = temperature the 
structure factor has a Dirac delta component at to = q = 
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S T =?{u J ,q)^(2^f5{u;)5{q) m 2 [6 U 



, <V,A 



(8) 



where the order parameter m is related to the expec- 
tation value of P^ L as (P' 1 ) = m^,^ — l/<5), the order 
parameter being aligned along direction jl. 

The lj = q = delta function is a characteristic of long 
range order and is absent for g > g c , where the structure 
factor has a delta peak at the quasiparticlc energy at 
T = temperature 

Sj=,V q) = A(g) (V - I/O) % - ?(?)) + .... (9) 

This behavior is schematically shown in Fig. [21 For Q < 
4 both the order parameter m(g) and the quasiparticlc 
residue A(g) scale to zero as we approach the quantum 
critical point g c , where the quasiparticle gap A vanishes 
and the quasiparticle description breaks down. 
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FIG. 2: (Color online) Schematics of the T = temperature 
structure factors on the paramagnetic (top) and ferromagnetic 
(bottom) sides of the phase transition. On the bottom we 
show the integrated structure factor, S(d) ~ J dq S(u>,q). 
Indices are dropped for simplicity. 

The T = temperature dynamics discussed so far 
gets drastically modified in the gapped phases at fi- 
nite temperatures where collisions between quasiparti- 
cles must also be taken into account. Fortunately, we 
can get an accurate description of the low temperature 
dynamics f or T <C A using the semiclassical formalism 
of Refs. Pl2rl3f| . 



In this limit, we can describe the spin system as a 
dilute gas of weakly interacting quasiparticles with mo- 
menta k ~ 0. Performing a perturbative calculation we 
find that, in the perturbative regime, for Q ^ 4 the 
scattering matrix S of these quasiparticles takes on a 
SU(Q — 1) -symmetrical form in agreement with the gen- 
eral arguments of Ref. ,26| . In the paramagnetic phase 
we find, e.g., 



6 A,A 



X' X X' 



(-1) 5* 6 



(10) 



thus quasiparticles exchange their internal quantum 
numbers under the collision while the many-body wave 
function picks up a phase (—1). 

Far on the ferromagnetic side, j < 1, we also find 
that the scattering matrix assumes the exchange form 
above, excepting that now can be viewed as kinks (do- 
main walls) between different vacuum states. In course of 
the scattering process these kinks exchange their quan- 
tum numbers in such a way that the the domain in the 
"middle" of the collision keeps its color (see Section HlBjl . 

Rather surprisingly, for Q = 3 the results found in the 
perturbative regime do not agree with the iS-matrix that 
one obtains by removing the cut-off and assuming that 
the system flows to a massive integrablc fixed point^i: 



(11) 



One, rather unlikely explanation is, that Eq. i|10|> only 
describes the scattering of q — * momentum quasipar- 
ticles far away from the quantum-critical point, and at 
for some finite values of g ^ g c 'phase transitions' occur 
where the asymptotic form of the quasiparticles jumps 
from one 'universal' form to the other. 

However, as we discuss in Appendix [CJ rather robust 
arguments constructed along the lines of Ref. |2(| sug- 
gest that Eq. Q10[l is indeed the (S-matrix of the q — > 
momentum quasiparticles on a lattice: The low energy ef- 
fective field theory contains three independent coupling 
constants, and unless two of these are the same, the <S- 
matrix takes the form (|l(Jfl rather than (|11|) . 

It is, nevertheless, quite concievable that the previ- 
ously mentioned two couplings become equal once one 
removes the cut-off A ~ j. 

This in turn means that for A < j an intermediate 
energy regime may exist, where, rather than Eq. (|10|) . the 
diagonal S- matrix Ijllfl of Ref. |3l] describes the scattering 
processes. 

In fact, numerical calculations of the mass spectrum for 
Q > 3 seem to support that the above Koberle-Swicca 
iS-matrix correctly describes quasiparticles at large mo- 
menta, q ~ A/c, where bound solitons are formed^ 

Therefore, while we believe that the correlation func- 
tions we derive based upon Eq. (|10|l always correctly de- 
scribe the regime T — > (q — > 0), the range of validity 
of the results in the regime \g — g c \ <C 1 needs further 
numerical investigation. 

We remark that there are a number of cases where 
in the 'universal scaling limit' important corrections are 
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missed that actually provide the leading contributions to 
some quantities, so that the Potts model would not be 
unique at all in this regard: 

A famous example is provided by the so-called free 
fermion point of the sine-Gordon model, where in the 
continuum limit one also finds a diagonal 5-matrix. 
There too, however, the asymptotic 5-matrix assumes 
the exchange form immediately once one introduces a 
cut-off or goes slightly away from this special point^ 
The same happens in the two-channel Kondo model in 
the scaling limit. There the channel susceptibility and 
also that associated with the composite superconducting 
order parameter simply vanishes in the scaling limit, be- 
cause the coefficients of these terms are inversely propor- 
tional to the cut-ofS^. There are also abundant examples 
in the literature where the physical properties of a model 
change dramatically once the integrability condition is 
violated (random matrix theory, integrable models with 
just a single impurity, etc.). 

The above simple form of the S- matrix, Eq. I|10(l . then 
enables us to compute analytically the correlation func- 
tion within the semiclassical approximation . This turns 
out to be a universal function of the typical separation 
£ c between the quasiparticles and their scattering time t, 
given explicitely by Eqs. 1(771)) and lf?TJ) . respectively. 

In the ferromagnetic phase we can express the spin- 
spin correlation function for T — > as 

(x, t) = m 2 1 (6^ - 1) R(x, t) , (12) 

where the relaxation function R is given by the follow- 
ing expression: 



JJ(x,t) = f ^e- |f| (1 - cos « (* e ~" 2+ucl (13) 



cos (sin(</>) x) 



(Q-i) 2 -i 



(Q-l) 2 + l + 2(Q-l)cos0 ' 



and x = x/£ c , t = t/r, and u = x/t, denote dimen- 
sionless separations, times and velocities. 

This formula is quite remarkable: On one hand, in the 
Q — > 2 limit it reproduces the exact results of Refs. 0, 
Il2j . However, while for Q = 2 the function R decays 
exponentially, for Q > 2 it has a diffusive structure for 
1 < t: 



R(x, t) ~ — = cxp ^ 
yt 



2t 



(14) 



This diffusive structure is related to the approximate 
SU(Q — 1) invariance of the scattering matrix. As a 
consequence, the T = temperature Dirac delta peak in 
the structure factor is replaced by 



From this, with the help of Q, we find that the suscep- 
tibility in the (semiclassical) limit, lo <^ T has a diffusion 
pole, 



1U 



QJ T Qq 2 - iur2^ 



(16) 



The prefactor in this expression, ~ m 2 £ c /T is just the 
static susceptibility, that can be interpreted as the Curie 
susceptibility of independent domains of size £ c having 
magnetization m. 

In fact, this diffusive structure is rather natural, and 
we shall discuss a simple explanation for it later. What is 
not natural, is the exponential decay found in the trans- 
verse field Ising model (Q = 2) J£ As we discuss later, this 
exponential decay is a consequence of long range correla- 
tions in the domain wall orientations along the chain for 
Q = 2, and is a peculiarity of the Ising model. 

For the low temperature correlation function on the 
paramagnetic side of the phase transition we find the 
following expression 



^i,ft' {.Xj t) 



A 5, 



1 



K(x,t) R(x,£) , (17) 



where the propagator K(x, t) is approximately the Feyn- 
man propagator of a quantum mechanical particle of 
mass A, and is given by 



K(x,t) 



-iA t 



A 

2-ni t 



exp 



A 



2c 2 t 



(18) 



for x and t within the light cone. Remarkably, the relax- 
ation function R{x, t) is the same as the one found in the 
ferromagnetic phase. This is very likely the consequence 
of the self-duality of the Q-state Potts modeliiSi 2 ^ 

We have to emphasize that Eqs. (|13f) and (|17f) rely 
on the structure (|ll)fl of the S-matrix. Therefore, while 
they are certainly valid in the regime far from the critical 
point, \g — g c | ~ 1, it is not clear below what tempera- 
ture these formulas describe the correlation functions for 
1 9 — <7 C | <C 1. This issue needs some further numerical 
investigation which is beyond the scope of the persent 
paper. 

It is important to emphasize that the results Eq. i|17|) 
obtained above are also relevant for S = 1 antifcrromag- 
netic spin chains; 2 ^^ and some two-leg ladder systems 
that can be mapped to a one-dimensional 0(3) quantum 
rotor model, also equivalent to the 0(3) sigma model in 
the long time and long wavelength limit i^li The one- 
dimensional quantum rotor model consists of a chain of 
ferromagnetically coupled quantum rotors, and is defined 
as 



5,,„> - — 



3„2 



(15) 



J 



ntn i+ i 



(19) 
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Here the vectors Hi denote iV-dimcnsional unit vectors 
and the Lj's denote the corresponding angular momen- 
tum operators. This model maps to the two-dimensional 
classical O(N) model. Correspondingly, it does not dis- 
play a quantum phase transition but has only a param- 
agnetic phase, where the coupling g always generates a 
finite gapA On the other hand, this paramagnetic phase 
is similar to that of the Q-state Potts model in that the 
quasiparticles have internal quantum numbers. Let us 
focus here on the probably most relevant N = 3 case, 
the 0(3) rotor model, where the gapped quasiparticles 
are spin S = 1 objects. Since the scattering matrix takes 
on the same universal form as in our case2&, all our cal- 
culations of the finite temperature properties carry over 
to this case as well and give 

(H{x, t) ■ H(0)) = A K(x, t) R Q=4 {x, t) , (20) 

where Rq = 4 is just the relaxation function Eq. (|13[) with 
Q = i. 

Thus our results are also relevant to spin S = 1 
antifcrromagnctic chains and also some of the experi- 
mentally studied spin ladder systems pi for which the 
one-dimensional 0(3) quantum rotor model provides a 
satisfactory description of low energy (long wavelength) 
fluctuations. 36 ! 37 

The scmiclassical formalism discussed above gives a 
consistent description of the dynamical fluctuations at 
very low temperatures in the gapped phases. Conceptu- 
ally, however, the most interesting regime is the quantum 
critical regime, A < T < j. In this regime dynamical 
correlations are governed by fluctuations related to the 
critical point, g = g c , and can be accessed by making use 
of the conformal invariancc of the critical theory. Confor- 
mal theory implies that the finite temperature dynamical 
susceptibility is given approximately by 

TV J • ttJ+eg N v rl ; u-cq \ 

x 1 V 15 _ 4ttT / 1 V 15 1 iirT > 

pf 14 _ Ai£±cq\ TV M _ 

A \15 1 AttT I L V15 1 iirT > 

for u>,cq,T > A, and correspondingly, the susceptibility 
exhibits u/T scaling in this quantum critical regime. A 
similar scaling form shall be obtained for the local sus- 
ceptibility 

Although we think they arc irrelevant for T — > 0, for 
completeness, let us also give here the relaxation func- 
tions obtained under the assumption of integrability in 
the gapped phases, i.e. using the diagonal S- matrix, 
Eq. In this case we obtain the expression 

RZz&i) = e-*^, (22) 
G(u) = -^e-" 2 + u erf (it) (23) 

for the decay function in the paramagnetic case, while in 
the ferromagnetic case we find 



R { g™{x,i) = e-i iG ^ . (24) 

In both cases the correlation function decays expo- 
nentially. The derivation of these expressions, which 
may be relevant in an intermediate temperature range, 
A > T > T* for \g — g c \ < 1, is given in Appendix |S1 

The rest of this paper is organized as follows: First, 
in Sec. [H] we investigate shortly the T = temperature 
properties of the system in the gapped phases by means of 
perturbation theory. In Sec. IIIII we use the semiclassical 
approximation to obtain the low temperature correlation 
functions. Sec. IIVI is devoted to the discussion of the 
quantum critical regime, and our final conclusions are 
summarized in Sec.0 Some of the technical details have 
been relegated to Appendices. 

II. PERTURB ATIVE ANALYSIS OF T = 
TEMPERATURE PROPERTIES 

In this section we shall study the scattering proper- 
ties of quasiparticles in the j > 1 and g -C 1 regimes, 
where straightforward perturbation theory allows one to 
obtain the energy of the quasiparticles, understand their 
structure, and determine their scattering matrix. 

A. Paramagnetic side, j > 1 

We shall first study the limit of very large g. In this 
limit we can perform an expansion in 1/g by taking the 
term describing the effect of transverse field ~ g as an 
unperturbed Hamiltonian, 

i 

and considering the ferromagnetic interaction as a per- 
turbation, 

i n 

The ground state of Eq. (|25|) is simply the paramagnetic 
state 

iAo)=ni A °>- ( 2? ) 

i 

where |Ao)j denotes the state |Ao)i <-> 1, . . . , 1) at 

site i in the basis of the ferromagnetic eigenvectors. El- 
ementary excitations in the g — ► oo limit consist of local 
'spin flips' of energy SE = jg, where one of the spins is in 
a P = eigenstate |A)», orthogonal to |Ao)j. Obviously, 
there are Q— 1 such states, and correspondingly, the label 
A runs from 1 to Q— 1. We shall denote the state having 
Potts spins flipped at sites ... by \i, A^; j, Xj; . . .). 
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Due to the ferromagnetic term, Eq. (|26fl . these local 
excitations can hop between lattice sites and get a dis- 
persion. In leading order in 1 / g the corresponding wave 
function and quasiparticle energy of an elementary exci- 
tation is given by 

A) « £>^ fc |j,A) , (28) 

3 

e\k) = e(fc)«j 5 (l- Alcosfc), (29) 
Q 9 

where k is the quasi- momentum of the excitations, Xj 
denotes the position of lattice site j, and we set the lattice 
constant to a = 1. 

It is also easy to compute the spin-spin correlation 
function in this approximation: The operator creates 
local spin flip excitations which then propagate along the 
chain. A simple calculation then gives a single-particle 
contribution to the dynamical susceptibility 

X* (Q, ^)t=o = (V -q) {LJ + i§) 2_ e 2 +-- > (3°) 

corresponding to a single particle contribution to the 
structure function given by Eq. @ with 

A(g)A~ A 4 / 15 ~( ff - 5c ) 2 / 9 (31) 

as one approaches the quantum critical point. 

Note that at T — temperature these quasiparticles 
have an infinite lifetime to all orders in the perturbation 
theory (as guaranteed by energy and momentum conser- 
vation and the existence of a gap). Of course, Eq. © only 
gives the leading behavior of S^' (u, q), and at higher 
energies a many-particle continuum also appears above a 
threshold due to higher order corrections to the ground 
state wave function. 

For j> 1 the perturbation theory in 1/g is convergent. 
Therefore the above picture holds to any order in 1/g for 
g > g Cl apart from the dispersion e(k) being renormalized 
and the quasiparticle weight A{g) being reduced for finite 
values of g. In general the single particle contribution to 
the correlation function reads, 

S^frt) = (<W - / £ A k ««*»-«<*>*> , (32) 

with the weight Ak related to the matrix elements of 
P M between the exact quasiparticle state |fc;A) and the 
ground state |0) as A k = ^ J2 x J{k; A 1 7^ =0 10) | 2 . 

For large values of x and t the integral in Eq. i|32|> 
can be evaluated within the saddle point approximation, 
yielding 

S£t' (x, t) = A(g) - K(x, t) , (33) 

with K(x,t) the Feynman propagator Eq. I|18|l and 
A(g) = linifc^o A k . 



The spin-spin correlation function further simplifies in 
the vicinity of the quantum critical point if we assume 
that the dispersion e(k) takes on a relativistically invari- 
ant form, 

e « \/A 2 + c 2 fc 2 , (34) 

with c being the quasiparticle velocity, and neglect the 
fc-dependence of A k — * A(g). In this case the zero tem- 
perature propagator can be expressed as 

K(x, t) = A K (A Vx 2 -cH 2 /c^ , (35) 

where Kq denotes the modified Bcssel function. For Q < 
4 the quasiparticle gap A and the weight A tend to zero 
continuously as one approaches the critical values of the 
coupling, g c . For Q > 4, on the other hand, the phase 
transition remains first order, and therefore both A{g) 
and A remain finite at the transition point. 

Having analyzed the single particle properties, let us 
now turn our attention to two-particle states. Wc look for 
a two-particle eigenstate of the Hamiltonian by making 
the following ansatz in leading order in 1/g: 

\k,k') = ^{vE^ I,+fc,Il) M;j,A')) 

A, A' i<j 

+ B XtX , ]T (e^+^li, A; j, A')) } . (36) 

i>j 

The coefficients A\ \i and B x ,\> in this ansatz are re- 
lated by the two-particle 5-matrix: 

A, A' 

Substituting l|H6^l into the Schrodingcr equation we ob- 
tain an equation for S x ' x , (k, k') that can be solved (see 

Appendix [X] for details). The expression of S xx , (k, k') 
is rather involved for general values of k and k' . How- 
ever, rather remarkably, for T — » 0, where states with 
both k and k' going to zero are relevant, S x x , (fc, k') re- 
duces to the simple form already given in the Introduc- 
tion, Eq. 10 : 

Thus the scattering matrix assumes the very same form 
as in Ref. j2||. While we obtained this equation in first 
order perturbation theory, we believe that it is valid to 
all orders in perturbation theory. 

B. Ferromagnetic side, j<1 

For j < 1 we can treat the ferromagnetic part of the 
Hamiltonian (|26f) as the unperturbed Hamiltonian, and 
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consider the 'transverse field' g, Eq. I|25|) as a perturba- 
tion. The ground state of Eq. I|26[) is Q-fold degenerate 
and corresponds to all spins being aligned in one of the 
Q possible directions, 

|0> M = J] Ms ■ (39) 

3 

For g — the excitations of Eq. 1261) consist of domain 
walls, 

iM,A=ni^n^v, (40) 

j<i i<j> 

and have energy j. The local field, Eq. (|25[) . generates a 
coherent motion of these domain walls. In leading order 
in g the wave function of the elementary excitations can 
be written as 

|fcW«5> to W)., (41) 

i 

and their energy is obtained by straightforward pertur- 
bation theory as 

#»'=j (l-g^cosk + ..)j . (42) 

Thus quasiparticlcs have a gap A(g) ks j(l — 2g/Q) in 
this phase too, and the ground state is stable. 

In the ground states the expectation value of the op- 
erators P M is finite, 

Correspondingly, the structure factor has a delta peak 
associated with the long range order at u> = q = 0, given 
by Eq. (jHJ). We remark here that any finite temperature 
induces a dilute gas of domain walls which destroy this 
long range order and broaden the delta peak. As we 
show in the next section, this delta peak has a Gaussian 
broadening for Q = 2, however, becomes a diffusive pole 
for any Q > 2. 

Similar to the paramagnetic phase discussed in the pre- 
vious subsection, all the above perturbative results carry 
over to any finite g < g c , because perturbation theory is 
convergent. For Q < 4 the order parameter vanishes as 
one approaches g c , 

rn(9)~(9c-9) P , (44) 

with (3 the order parameter exponent for the two- 
dimensional classical Q = 3 state Potts model (fe = 
1/9) and the two-dimensional Ising model (f3 2 = 1/8), 
respectively^ and the quasiparticle gap vanishes with 
the same exponent as in the paramagnetic phase, v = 5/6 
for Q = 3 and v = 1 for Q = 2. 

We now turn to the study two-particle properties in the 
limit g <C 1. To keep track of domain wall excitations, 




FIG. 3: (Color online) Sketch of colliding domain walls on 
the ferromagnetic side at T <C A as k — > 0. Note that the 
polarization of the middle domain must be conserved and so 
the order of the domains along the chain does not change with 
time. 



it is worth labelling multi-domain wall configurations in 
a slightly different way as before. Suppose that we have 
a domain wall with a polarization jj, on its left and yj 
on the right. We can then define the quantum number 
6 = 1, . . . , Q — 1 of this domain wall as the size of the 
step between the two sides, 

fx' = (/i + 6>) mod Q . (45) 

Clearly, to characterize any configuration, it is sufficient 
to give the vacuum state on the left of the chain, and then 
specify the quantum numbers {6j}. With this notation, 
we can thus denote the state described by Eq. I|4U|) as 

|M,A = N,6V (46) 

With the above notation, the two-particle wave func- 
tion has the following form in leading order in g, 

\k,k') = Y J { A ^ , Y,{ ei{kXi+k ' Xi) \ i ^^ e ')') + 

6,8' i<j 

" '•":./•"' ,)}• (47) 

i>j 

Again, a simple calculation outlined in Appendix ^ pro- 
vides us the scattering matrix relating the ampli- 
tudes Agfii and Bg gi of the incoming and outgoing par- 
ticles, respectively. In the limit of vanishing momenta, 
k, k' — > we find 

i.e. quasiparticles scatter as 'hard balls'. This equation 
can also be visualized condition that the orienta- 
tion of the order parameter between two colliding domain 
walls with vanishing quasi-momenta remains unchanged 
after the collision. Again, we believe that the exchange 
form of the scattering matrix, Eq. (|48|l . remains valid to 
all orders in g for g < g c . 
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III. SEMICLASSICAL APPROXIMATION IN 
THE GAPPED PHASES 

In this section, we shall study dynamical correlation 
functions in the T <C A limit in the gapp ed phases , 
following the semiclassical treatment of Refs. |!2l26l30l| . 
This approach is based on the observation that at very 
low temperatures only quasiparticlcs with energy close to 
the quasiparticle gap are present, e ~ A. The energy of 
these quasiparticlcs can be approximated as 

£fc = A+ A~y + ° (fc4) ' (49) 

with c a constant playing the role of the speed of light, 
and A/c 2 the mass of the quasiparticlcs. For T <C A 
the distribution of quasiparticlcs is described simply by 
Boltzmann statistics, 

n(k) ~ e-" A e-P ^ , (50) 

with (3 = l/T the inverse temperature. Correspondingly, 
the quasiparticle density is exponentially small 

•-w- 1 '^'" 4 "- <51) 

and the typical separation between them, dx ~ 1/ ' Q in- 
creases exponentially at low temperatures, dr ~ e A / T . 
This must be compared to the De Broglie wavelength of 
the particles, At, giving the quantummechanical exten- 
sion of the quasiparticles' wave function. This latter is 
given by the inverse of the typical momentum k of the 
quasiparticlcs, At ~ c/VtA. Clearly, at very low tem- 
peratures the average separation of the quasiparticles is 
much bigger than their quantum mechanical size, 

d T > At , (52) 

which makes one possible to treat the quantum mechan- 
ical state of the system in both phases within the semi- 
classical approximation for T <C A. 

Unfortunately, in one dimension neighboring particles 
cannot avoid collisions with each other, thus they will 
get unavoidably closer than the De Broglie wavelength, 
where quantum mechanics is at work. These collisions 
must therefore be described within the framework of 
quantum mechanics. Fortunately, at low T, where the 
semiclassical limit is valid, the system is dilute enough so 
that we have only 2-particle scattering. Furthermore, the 
colliding quasiparticles have momenta k ~ 1/c?t- There- 
fore, at low enough temperatures, we can use the simple 
forms of the scattering matrices given by Eqs. I|38[) and 
1)48(1. This simple form of the 5-matrix will enable us 
to obtain analytical results for the correlation functions 
within the semiclassical picture. 



A. Semiclassical correlation function on the 
ferromagnetic side 

Let us first compute the spin-spin correlation in the less 
complicated case of ferromagnetic phase. By definition, 
the time-dependent correlation function is defined as 

S^,(x,t) = (e itff P"(a;) e~ itH P>'(0)) , (53) 

where (...) denotes now thermal averaging over all pos- 
sible many-body states, 

(...)=5>|..»e-^», (54) 

n 

where the summation runs over all many-body eigen- 
states \n) of energy E n of the total Hamiltonian. Within 
the semiclassical approximation we can replace this aver- 
age by an average over all possible quasiparticle configu- 
rations, i.e., by an average over the quasiparticle veloci- 
ties v Vl their positions x Vl and internal quantum numbers 
V , (y = 1, . . . , M with M the number of quasiparticles), 

(...) « ]T fl[dx v l[dv u p({x u ,v v ,e u }) 

{M v v 

{{x v , v v , 6 V }\... \{x v , v v , 9 V }) , (55) 
where the distribution function P({x u , v v , #„}) is simply 

P({xu,v u ,e v }) = ^f ll%)- ( 56 ) 

with L the system size and -P(vp) the Boltzmann distri- 
bution of the quasiparticle velocities v = c 2 k/ A, 

P(v) = 6 " ^ " • (5?) 

Of course, in the above equations an average over M 
should be taken. In the L — * oo limit, however, we can 
replace M by the average particle number, M —> g L, 
without changing the final result. 




FIG. 4: (Color online) Time evolution of domain walls on 
the ferromagnetic side at T <C A. For the correlation func- 
tion S^^i (x, i) we need to take the average over all domain 
configurations. 

To evaluate the correlation function S M(I ^s,t), let us 
first compute the probability Qnn>(x,t), that at time 
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t = and position x = the order parameter points 
into direction /j while at time t and position x it points 
into direction To evaluate this probability, we notice 
that the number of domain walls is conserved in course 
of every domain wall collision. Furthermore, the 'color' 
(orientation) of the g'th domain from the left remains 
unchanged due to the exchange form of the scattering 
matrix (see Fig. 0}. As a consequence, Q^^{x,t) can be 
simply written as 



Q„ tl >>(x,t) = 5^P((0,0) in Z'th) 



Lk 



P((x,t) in fc'th) P(l,n\k,n') 



(58) 



i.e. it is simply a product of the probability 
P((0,0) in Z'th) that (t = 0,x = 0) is within the Z'th 
domain from the left, while the (i, x) point lies in the 
fc'th domain, and probability that the Z'th domain points 
in direction « while the fc'th in direction fjf, P(l, fi; fc, fjf). 
We have to remark here that the orientation of domain 
walls is correlated. This correlation is simply generated 
by the constraint that if a domain points in a direction « 
then its neighbor must point into one of the other Q—l di- 
rections. This correlation is hidden in P(l, u; k, fjf), which 
can be computed by constructing a Markov equation, and 
depends only on n = I — k, P(l, u; k, fjf) = P{n, fi, fjf): 



P(n,fJt,fjf) 



1 



1 



Q \Q- 1 



M 



1 



(59) 



Note that this correlation function decays exponentially 
for any Q > 2, while it oscillates for the Ising model, 
Q = 2. 

The first part of Eq. (|58|l also depends only on the sep- 
aration n = l — k and the coordinates x and t. Therefore, 
Eq. (511) can also be rewritten as 



(x,t) 



oo 



D(n, (x,t)) P(n,n,ijf) 



(60) 



with P(n, fi, ft') given by Eq. (|59|) and 

D(n, (x, <)) = p ((°> °) in l'th)P((x, t) in I + n'th) . 

(61) 

Clearly, D(n, (x, t)) is just the probability that the do- 
main of (x, t) is the n'th domain to the right from the 
domain of (0,0). 

The probability D(n, (x, t)) can be computed as fol- 
lows. First, following Rcfs. let us introduce the 
notion of 'particles'. In the configuration space, domain 
walls are located along straight lines v, 



x v (t) = x u + iv t , 



(62) 



where x v and v u are the position and velocity of the i/th 
domain wall at t = 0. 'Particles', however, correspond 
to a given step 9 and arc 'reflected' when two lines cross. 



There is a simple way to tell which particle p is moving 
along line v at time t. This is given by the function p v {t) 
which we chose to coincide at t = with the index v of a 
line, p u {t = 0) = v. Particles are, however, impenetrable. 
Therefore if line v is crossed by another line from the left 
then p v decreases by one, while if it is crossed by a line 
from the right, then it increases by one. Therefore, to 
keep track of p v (t) we just have to count the number of 
lines that crossed x v (t) from the left and from the right 
since time t = 0, 

Pu {t) = [©(**(*) - M*)) - - av)] . ( 63 ) 

v' 

where denotes the step function. 

Similarly, for a fixed set of lines, {x u , v u }, the probabil- 
ity that the straight line from (0, 0) to (x, t) moves across 
?i domains to the right is simply 

P{n\{x u , v v ], (x, t)) = S n j2je(x-x v (t))-e(- Xv )] > ( 64 ) 

To obtain D{n 1 (x,t)), we should average this expression 
over {x„,tv}. To do this we just introduce the following 
integral representation of the Kroneckcr-dclta, 



*n,^ K [e(x-x„(t))-e(-x„)] 



# M (n-V [&(x- X „(t))-e(-x„)]) 

2tt 



(65) 



The advantage of this form is that the averaging fac- 
torizes and can be evaluated analytically. With some 
algebra we obtain that 



D(n, (x,t)) 



2tt 



e m,t >I(x, t) 



I(x,t) = (e-WeCs-^-^tJ-ec-*,,)}) 



(66) 
(67) 



The integral I(x, t) can then be simply evaluated to yield 

I(x, t) = 1--L {G(u) + G(-u) - e l(t >G{u) - e^G(-u)} 

(68) 



where 



G(u) 



2V^ 



crfc(w); 



(69) 



and we introduced dimcnsionless time and length, x, = 
x/£ c , and i = t/r and the corresponding dimcnsionless 
velocity u = x/i, with £ c and r the characteristic classical 
correlation length and time 



£ c = i/q 



i 



Q - 1 V TA 



27TC 2 



A/T 



1 



A/T 



(70) 

Finally, I(x, t) M in Eq. (JHEJ can be re-exponentiated as 

I(x,t) M = e -t{G(u)+G(-u)-e^G(u)-e-"*-G(-u)}^ /^n 
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Substituting Eqs. (JBSJ and into Eq. lO the sum 

over n can be computed using the identity 



E 



i 



»i 



(Q-i) 2 -i 



(Q-l) 2 + l + 2(Q-l)cos0 
(73) 

and the final form of the function Q^^'{x, t) reads: 



Q- 



(74) 



with the relaxation function given by Eq. (|13|) . The cor- 
relation function S^^i{x,t) can then be easily computed 
in terms of Q^^'{x^t) using Eq. i|43|) . as 



A, A' 

Putting all this together we finally obtain 

/ 1 \ 



(75) 



(76) 



0.05 



0.04 



0.03- 



0.01 




0.02 



200 



FIG. 5: (Color online) The relaxation function R(x, t) as a 
function of t for different values of x at Q = 3. 

The time evolution of R(x,t) is shown in Fig(5]The 
behavior of the relaxation function R(x, t) is quite inter- 
esting for t ^ 1. The integrand in the definition 113|) 
goes to zero exponentially as t — > oo, and only a small 
range around cf> f» contributes to the integral. Hence 
the relaxation function can be approximated as 



R(x, t > 1) 



d<P 
2tt Q 



(Q-2) 



cos(o>a:) , 
(77) 

where F(u) = l/y^e - " + u erf(u). 

This integral can be calculated analytically and is given 

by 



R{x,t^> 1) 



Q 



i 



(78) 



For f > i we have -F(it) ~ l/v^, an d the relaxation 
function assumes a diffusive form 



i?(a; < t) cx 



Co 



-a; 2 /4Dt 



(79) 



with the diffusion constant defined as 

i e c 



D 



2V* r 



(80) 



As a result, the dynamical structure factor has a diffu- 
sive structure at u),q — > 0, as already discussed in the 
Introduction. 



B. Semiclassical dynamics in the quantum 
paramagnetic phase 

The analysis of the previous subsection can also be 
extended to the paramagnetic phase. Similar to the 
ferromagnetic phase, the average over the many-body 
eigenstates in Eq. (|54|) can be replaced in the semiclas- 
sical limit by an average over all possible initial states, 
\{x v , v v , ■ The main difference in the calculation 
is that on the paramagnetic side the operator P^(0) in 
Eq. H53f) creates a quasiparticle with some internal quan- 
tum number Ao, and velocity v at time t — with some 
probability amplitude e^ g (y) in addition to the already 
existing quasiparticlcs. This particle, together with the 
other quasiparticlcs, propagates under the action of H 
in Eq. (|53|) ) and collides with them. For very small 
temperatures, however, the two-particle scattering ma- 
trix takes on an exchange form, and therefore particles 
only exchange their velocity while conserving their inter- 
nal quantum numbers. 



t = 



t = 



backward 
evolution 



forward 
evolution 



x= 



x 



FIG. 6: (Color online) To evaluate the correlation function 
on the paramagnetic side, one must keep track of the parti- 
cle created by the operator . Lines drown with different 
line-styles denote trajectories of different 'particles'. The two 
states after and before the forward and backward time evo- 
lution must be identical. Therefore, the internal quantum 
numbers of the particles must obey constraints. 

As a consequence, at any time precisely one of the par- 
ticles will have the velocity v of the initial particle and it 
will be at position vt. It is easy to see that this very parti- 
cle must be annihilated at time t by P^ (x) otherwise the 
final state obtained after the action of e ztH in Eq. I|53|) will 
be orthogonal to the initial state. (The evolution of these 
quasiparticles is shown in Fig. 0) The probability am- 
plitude that this particle is annihilated is proportional to 
( e i'W)* e ' fa with k = A v/c 2 the quasi- momentum that 
corresponds to v and A the internal quantum number of 
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the particle that is removed. As we shall see below, A' in 
this expression must be equal to A . This shall be guar- 
anteed by another requirement, namely that the internal 
quantum numbers in the final state must be exactly the 
same as those of the initial state. To see what this con- 
dition implies let us assume that in the initial state we 
had quasiparticles with quantum numbers {Ai, A2, • ■ •} at 
positions x^ [x\ < a; 2 < . . .) with velocities v^, and that 
the new particle is created right after the p'th particle 

P"(0): {...,A P ,A P+1 ,...} (81) 
1 — ► {. . . , X p , Ao, Ap+i, . . .} . 

Let us now define the function po(t) which gives us the 
label of that particle which has velocity v at time t. Ob- 
viously, po(t = 0) = 0. The important observation is that 
the above order of these quantum numbers from the left 
to the right along the chain does not change under the 
collisions due to the simple form of the scattering matrix 
Eq. I|38|) . Suppose now that after the time evolution t 
particle 'n + p 1 moves with velocity v, i.e., Po[t) takes 
the value po(t) = p + n. (Here we assume for the sake 
of simplicity that n > 0, but the derivation can be eas- 
ily generalized for n < 0.) Since the operator P M must 
destroy precisely this particle of velocity v, P M therefore 
changes the series of quantum numbers as 

P M : {. . . , X p , Ao, Ap+i, . . . (82) 

. . . , Ap+ n — 1, Xp-\-m Xp-\-n-\-l: ■ • ■} 

1 ► {. . . , A p , Ao, Xp+i, ■ ■ ■ 
. . . , Xp-\~ n —i , Ap_|_7 t _|_ 1 , . . .} . 

Since the order of these quantum numbers does not 
change under the action of e~ ltH either, requiring that 
the internal quantum numbers of the initial states be the 



S^(x,t) ={j2jdv «(«))*<'(«) e -*(*Mx-*0«) 



where the factor 1/(Q — 1)'™' comes from the condition 
Eq. H83|) after averaging over all possible internal quan- 
tum numbers. The first term in this expression gives just 
the T — correlation function Eq. i|33|) , while the second 
part is the same relaxation function as the one found in 
the ferromagnetic phase, Eq. (|13fl . 

S^,(x,t) = S^ (x,t) R(x,t) . (87) 

Finally, let us briefly discuss the properties of quasi- 
particle density-density correlations in the paramagnetic 
phase. The density of quasiparticles with quantum num- 



same as those of the final state amounts in the condition 

■^0 = ^p+1 = • • • = ^p+n ■ (83) 

Now, similar to the ferromagnetic case, the value of n in 
the previous expressions can be simply computed as the 
number N + of lines x v (t) that cross the line connecting 
(0, 0) with [x, t) from the right minus the number of lines 
A_ that cross it from the left, 

n = 5}e(x(t) - x«(t)) - 6(0 - x v )] , (84) 

where v labels the quasiparticles in the initial state, 
x v {t) = x v + v v t, and x(t) = vt is the trajectory of the 
new particle with velocity v is created by P M (0) at the 
origin, x(Q) = 0. 

To compute the contribution of this particular state to 
the correlation function we must consider the important 
detail of phase factors. Quasiparticles in the initial state 

generate a phase factor e~ lt ^^ e '""' ) under the action of 
e -itH _ This is, however, completely canceled under the 
action of e ltH , except for the quasiparticle created by 
P^ giving a factor e -' lte ( v ) . Furthermore, every collision 
results in a sign change of the many-body wave function. 
All these signs cancel under the forward and backward 
propagation, excepting the ones that are associated with 
collisions with the extra particle. These give an extra sign 
(— 1) which can, however be also written more 

conveniently as 

S= (-i) N +- N - = (-1)™ . (85) 

Putting all these together we obtain the following ex- 
pression for the correlation function, 



iq_ i)|ni a ",£je(*-M*))-e(o-z„)]) { , ' ( 86 ) 

— — OO ^ ' V ' U 

I 

ber A is simply given by 

Q X (x,t) = J2Sx.x P:At) 5(x u (t)~x) , (88) 

where the summation runs over all lines. The average 
density, g = (g(x,t)) is simply given by Eq. (|STJl . The 
density-density correlation function can be computed an- 
alytically within the semiclassical approximation. We 
shall not discuss this calculation here, since it is basically 
identical to the computation of Damle and SachdevS 
The result is simply 
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(g\(x,t)gx>(0)) = g 2 



OA, A' 

Q-l 



e -t (G + +G_) x 




+ 2C+C\ I (2ty/G^Gl) + 



G\ G- + C 2 _ G H 

yj G+G- 



h{2t^G^GZ)^ 
(89) 



where G± and C± denote the functions 

G ± = G(±u) = J- e"" 2 T ~ erfc(±u) , (90) 
2V7T 2 



C± = C(±w) = - erfc(±ti) , 



(91) 



with u = x/i the dimensionless velocity defined earlier. 

For large time scales, t 3> 1 the density-density corre- 
lation function also displays a diffusive behavior^ 



(g\(x,t)g y (0)) ~ oa, A' e" 



-0?s 2 /2* 2 i (92) 



in agreement with the behavior of the relaxation func- 
tion R(x,t). In the Z% case, this correlation function 
also bears physical meaning and corresponds to the cor- 
relation function of the chirality density. 



IV. DYNAMICS IN THE QUANTUM CRITICAL 
REGION 

In this section we shall focus to the Q — 3 Potts 
model which has a quantum critical point at T = and 
.9 = 9c- We can gain many information already from 
the quantum-classical mapping, by just using the scaling 
properties of the singular part of the free energy density: 

f(g - g c , h, T) = b- 2 /(6*«(5 - g c ),V*h, Tb) , (93) 

where yt = 2 — Xt and y^ = 2 — Xh denote the scaling di- 
mensions of the temperature and the magnetic field in the 
classical Potts model, and x t = 4/5 and Xh = 2/15 de- 
note the dimensions of the corresponding primary fields 
that are known from conformal field theory>i2i2i Note 
that the temperature plays the role of a finite system size 
while the coupling g corresponds to the temperature. 
From Eq. I|93|l immediately follows that the gap van- 



ishes as A 



1.9 



15/6 



, and that the susceptibility 



behaves in the quantum critical region as 



X(T > A) 



1 



^26/15 



(94) 



For T < A the susceptibility becomes finite, but it 
diverges as one approaches the quantum critical point, 



X(T = 0) 



1 



I. 9 -.9c | 13/9 



(95) 



So far, we only discussed thermodynamical properties. 
However, a lot of information can also obtained by mak- 
ing use of the conformal invariance of the critical theory. 



At the critical point, the imaginary time correlation func- 
tion is scale invariant and has the following form: 

x£^W) = (S^-hx(x,r), (96) 



1 



(t 2 + X 2 ) Xh 



(97) 



For simplicity, we set the 'speed of light' to one in this 
section. 

To obtain the finite temperature imaginary time sus- 
ceptibility (retarded correlation function), we introduce 
new, complex coordinates z = r + i x and z = r — i x 
and rewrite l|§7|) as 

x r=°( Zl z)~_i_ . (98) 

z z 

Then we map the complex plane to a strip of width f3 = 
1 /T by the transformation 



w = cotan (n T z) 

it T 



(99) 



and use the transformation properties of primary fields 
to obtain^ 



x(x,t) 



(sin(7rT(r — i x))) Xh (sin(7rT(r + i x))) Xh 

(100) 

To obtain the (retarded) susceptibility x( w ) l) one nas to 
Fourier transform this correlation function to obtain the 
Matsubara Green's function, x(w„,g) and then analyti- 
cally continue it back to the real axisi^ The suscepti- 
bility has the same structure as the one obtained for the 
transverse field Ising modelji 



1 



L V 15 



' AttT ' 



rf — - i— 

1 V 15 1 4tt 



T ' 



7^26/15 TV 14 _ ■ u+q \ TV 14 _ - Ul-q \ 
1 V 15 AttT > L \15 1 AttT' 

This immediately implies that the dynamical suscepti- 
bility x( w ) = xi^iQ = 0) shows ui/T scaling, 

C 



F(u/T) , 



y26/15 

where the scaling function F(y) is given by 

F(y) 



r(&)r(8-<£), 

and has the following asymptotic properties: 
f 1 



r(i4/i5) 
r(i/i5) 



-26/15 i 7t 13/15 



(101) 



(102) 



y < i , 

(103) 
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Another quantity of interest is the local susceptibility, 
X loc (w) = x{ u i x = 0)i tnat can a l so be measured by 
neutron scattering, and which can be computed to behave 
as 



C 



X l ° cal (co,T) = lflm G(uj/T) 
with the scaling function 



G(y) 



r(4)r(g 



(104) 



(105) 



'lit' 



having the asymptotic properties: 



G(y) " { Wiw (y/ 2 *r 11/15 ^ * 11/30 v >> 1 • 

(106) 

We have to emphasize that, although conformal invari- 
ance only holds at the quantum critical point, the above 
expressions also apply to the entire quantum critical re- 
gion, j 3> T,uj,q ^> A, where the quasiparticle gap does 
not play an important role. 



V. CONCLUSIONS 

In this paper we studied dynamical properties in the 
gapped phases of the Q-state Potts model in the T — > 
limit, and in its quantum critical behavior for Q = 3 
(the Q = 2 Ising case has already been studied thor- 
oughly in the literature). Deep in the gapped phases the 
T — > correlation functions were found to show a diffu- 
sive character for Q > 2. This is a consequence of the 
simple structure of the low energy scattering matrices, 
, and can be understood in a very simple way. Con- 
sider, for the sake of simplicity, the ferromagnetic side, 
T«A,j< g c . In this limit domain walls propagate as 
particles and eventually collide with other domain walls. 
However, due to the simple structure of the scattering 
matrix the order of the color of domains does not change 
under these collisions. 

Therefore, if we look at a given domain, then its size 
and position changes in time, but its color (orientation) 
does not. In other words, each domain has a typical 
size ~ £ c and a given domain diffuses along the chain as 
time flows. The colors of domains far away from each 
other are uncorrelated. Therefore, the probability that 
the domain at time t at position x has the same color as 
the domain at t = and x = is approximately given by 
the probability that the domain wall that was at position 
x = at time t — diffused to position x under time t, 
and is therefore proportional to 



1 



,-x 2 /4,Dt 



(107) 



with D the diffusion constant of domains, given analyti- 
cally by Eq. (glty. 



This simple argument, however, obviously fails for 
Q = 2, where the colors of domains far away from each- 
other are strongly correlated, as illustrated in Fig. 
Mathematically, from Eq. I|59|l we see that correlations 
between the colors of two domains separated by n domain 
walls decay exponentially for Q > 2, C ~ e~™ MQ- 1 )^ 
while they simply oscillate for Q = 2. These correlations 
for Q — 2 lead to a destructive 'interference' and give 
rise to an exponential decay of the correlation functions 
for the Ising case, Q = 2ml Note that the Q = 2 Potts 
model is integrable^, and a similar exponential decay 
has also been found recently by Altshuler and Tsvelik in 
some other exactly integrable cases ml 

On the other hand, after this work was completed, 
we became aware of independent work by Damle and 
Sachdev^S who find that diffusive behavior appears in 
the correlators of the one-dimensional sine-Gordon model 
too. Rather remarkably, the decay functions in the two 
cases seem to differ only in the choice of the parameter 
Q — 1, which takes a value Q — 1 — > — 1/ cos(27r?7/7) in 
the sine-Gordon calculation, although the two models are 
completely different. This similarity between the two re- 
sults indicates that the form of correlation function we 
have found is remarkably universal. 

We have to emphasize that the diffusive character of 
the correlation functions for Q > 2 is a consequence of 
the simple structure l|10|) of the scattering matrix be- 
tween quasiparticles of vanishingly small momenta. At 
finite temperature, however, the colliding particles have 
finite momenta ~ yTA/?, and therefore the quantum 
numbers of the colliding particles change with a finite 
probability in course of a collision. As a consequence, 
the color of the domain between the two colliding do- 
main walls can change with a finite probability. At long 
enough time scales, this process should lead to a decay of 
correlations, unless the correlation function of the density 
of a conserved quantity is computed^ 

The probability of color changing can be estimated to 
be proportional to ~ (a 2 (Ap) 2 ) with Ap the typical mo- 
mentum difference of the two incoming particles and a 
the lattice constant, and scales as Pfu p ~ a 2 A T/c 2 . A 
given domain shrinks to zero due to the motion of do- 
main walls approximately t/r = i times during a time 
period t, and each time it shrinks to zero it changes color 
with the above probability. 

As a result, the diffusive form of the decay function 
in the quantum Potts model must break down at a time 
scale 



idiff 



a 2 A T 



(108) 



above which correlation functions must decay exponen- 
tially. 

We also present very robust renormalization group ar- 
guments to show that on the lattice the asymptotic <S- 
matrix gencrically takes the exchange form found in the 
perturbative regimes, and thus the diffusive correlation 
function we found should also describe the vicinity of 
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Ising model 




Potts model 

12 3 1 2 13 

FIG. 7: Domain structure in the transverse field Ising model 
and in the Q — 3 state quantum Potts model. For Q = 2 
orientations of far away domains are correlated while for Q > 
2 they are not. 

the quantum critical point as T — > 0, \g — g c \ <C 1. In 
this regime, however, further numerical calculations are 
needed to determine the applicability and range of valid- 
ity of our formulas. 

Our results are not valid for the case Q = 4, where scat- 
tering matrix cannot be computed in the way outlined in 
Appendix [H] because the matrices become singular for 
k, k' — > 0, and cannot be inverted. As a consequence, the 
scattering matrix has a singular structure for k,k' — > 0, 
and a special treatment is needed. Indeed, the Q = 4 
state Potts model is known to have unusual thermody- 
namical properties due to the presence of a marginal op- 
erator at the critical point^ 

Also, in the present paper we neglected umklapp pro- 
cesses. It has been shown recently, that in one dimension 
such processes may dominate the relaxation of some cur- 
rents that overlap with other conserved quantities^ 

These processes are, however, probably not relevant for 
the correlation functions computed here, since "spin con- 
servation" in course of the scattering process is anyway 
only approximate in the Potts model. These umklapp 
processes are, however, probably important for energy 
relaxation, which has not been studied in this paper. 

Finally, let us discuss physical applications of our cal- 
culations. Though the 3-state quantum Potts model 
can be realized in trigonal ferromagnets in a magnetic 
field, and also recently some spin chain models have been 
shown to map on this model^ it is physically not easy 
to find realizations of the 1-dimcnsional Q-state quan- 
tum Potts model. However, our results obtained for 
the paramagnetic phase are very general, and carry over 
to gapped antiferromagnets and spin ladders with little 
modifications. For the spin-spin correlation function of a 
one S = 1 spin Heiscnbcrg antifcrromagnet, e.g., our re- 
sults imply that the spin-spin correlation function decays 
as 

(S( Xi ,t) ■ S(0)) = (-1) 4 AK{x t ,t)R Q=i {x t ,t) , (109) 

with the Feynman propagator K and the decay function 
R given by equations l|18fl and i|13|) . This result implies 
that the dynamical susceptibility has an approximate dif- 
fusive structure also at momentum q = it /a, similar to 



the diffusive structure that appears at q = 0i2& However, 
while the q = diffusion pole follows from the SU(2)- 
invariance of the Hamiltonian and is therefore protected 
by symmetry, the q = it /a diffusive character is only ap- 
proximate, and is a consequence of the simple exchange 
form of the scattering matrix. 

We would like to express our special thanks to Kedar 
Damlc and Subir Sachdev for showing us their preprint 
before making it public, and F. Igloi for valuable com- 
ments and discussions. We would also like to thank 
A. Tsvelik and R. Konik for useful discussions. This 
research has been supported by Hungarian Grants No. 
OTKA T046303, T046267, and NF061726. 



APPENDIX A: DERIVATION OF THE 
SCATTERING MATRICES 

1. Scattering matrices - ferromagnetic side 

To obtain the two-particle scattering matrix in leading 
order in g we make the simple ansatz Eq. (|47|l for the 
two-particle wave function |^), 

I*) = E { A w E (e* (te<+fc ' xj) IM; + 

6,6' i<j 

Be^^+^iAj^h) } . (Al) 

i>j 

where jl is the orientation of the chain on the far left, and 
9 and 9' denote the two kinks corresponding to the two 
domain walls. In leading order in g and for coordinates 
\i — j\ 1, H g just moves the two domain walls inde- 
pendently, and |<F) is clearly an eigenstate of the total 
Hamiltonian with an eigenvalue 

H\*) = (E + e(k) + e(k'))\V) , (A2) 

Eq being the ground state energy, and e(fc) the quasipar- 
ticle energy given by Eq. (|42|l . However, |\&) must satisfy 
Eq. I|A2() also for j = i + 1, i.e. for nearest neighbors. 
Observing that the operator H g = —jg J2 n Pn just flips 
each spin to some other direction, we can write in the 
original notation 

n 

(/i|P|/x) |/2, i- l,fj,,i + l,fi') + 
(li\P\l/) |/2, i,/M + 2,//) + 
]T (//'!%) \ fi,i,^",i + l,n')+-- - , (A3) 

where 9 = (fj, — /2) mod Q, 9' = (// — /t) mod Q, and we 
neglected all other terms involving more than two domain 
wall excitations, since these are high up in energy. Since 
the off-diagonal matrix elements of P n are all equal to 
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1/Q, we can write this in the 'kink' representation as: 
1 + 1,0%= (A4) 

^ [h'-l,0; i + l,0% + M; i + 2,0% + 



where j denotes the Kronecker delta modulo Q. Pro- 
jecting out the i + 1 = j component of the Schrodinger 
equation Eq. (|A2|) we obtain the following constraint for 
the coefficients Ag g> and Be g>, 



E 



s X + 



i 



gifc' _i_ g— ik 0+6 



E 



<5 2<s ; + 



This equation can clearly be inverted to give the two- 
particle 5-matrix in leading order in g, but the solution 
is rather complicated even for Q = 3. Eq. <| A5|> simplifies, 
however, in the limit fc, k' — > 0, relevant for very small 
temperatures, T <C A. 

For Q = 3 we obtain in this way 



s e e: e ,(k,k' ^o) = (-!)<$ 4 



(A5) 



We remark here that the above result holds for any 
Q 7^ 4. The Q = 4 case, however, seems to be special: 
then the operator in front of the coefficients Ag g, and 
B$fii in Eq. (|A5|) has zero eigenvalues for k = k' = 0, the 
inversion is problematic, and the 5-matrix does not take 
the form HASH . 



2. Scattering matrices - paramagnetic side 

To obtain the two-particle scattering matrix in leading 
order in 1/g we follow similar steps as for the ferromag- 
netic case. The ansatz for the two-particle wave function 
can be written as follows: 



I*) = E { A W E (e? (kxt+k ' Xi) \i, A; j, A') 

A, A' i<j 

Biv^^'^M;^'))}. (A6) 

For \i — j\ 3> 1 |$) must satisfy a similar two-particle 
Schrodinger equation as l|A2f) . and so it must be valid 
for j = i + 1 too. Let us now calculate the effect of 

iZferro = -j En Eu P n P n+l 0n tnC * + 1 = 3 term 01 

|\]/). If we neglect the high energy terms with multiple 
particles we get 



EEW+i M;i + 1,A') = 



(A7) 



1/Q li-l.Ajt + l.A') + 
1/Q |i,A;i + 2,A') + 

^^<A'|^|A)<A|^|A')N,A'; l + l,A). 



XX' M 



Substituting this to the two particle Schrodinger equa- 
tion leads to a similar constraint for the coefficients as in 
the ferromagnetic case: 



E 



A, A' 



i 



-r, t m xx' 



X'X 



(A8) 



E 

A, A' 



^A.A' ) 



where 



=^(A'|P"|A)(A|P^|A'}. 



(A9) 



Similar to the ferromagnetic case, in the limit of van- 
ishing quasiparticle momenta k,k' — ► 0, this equation can 
be solved to obtain the 5-matrix for Q = 3: 



A' rA 



(A10) 



APPENDIX B: COMPUTATION OF THE 
RELAXATION FUNCTIONS ASSUMING A 
DIAGONAL S-MATRIX FOR Q = 3 

1. case: Ferromagnetic side 

Let us compute the correlation function C f °"°(x, t) de- 
fined as the probability that the domain containing (0, 0) 
has the orientation fj. while that at (x, t) has orientation 
y! . Note that C w < is defined in terms of the P^'s rather 
than the P^s. 

We can calculate this probability by just keeping track 
of the number iV + and AT_ of kink excitations (domain 
walls) with quantum numbers + and — , respectively, that 
cross the line between the points (0,0) and (x,t). The 
points (0, 0) and (x, t) have the same domain orientation 
if N + equals N_ up to modulo three. 

Thus the correlation function is 

^T°(x,f) = i^P n [(0,0)^(x,t)]x 

n 

x E(fc)^( 2fc -™ + A ^)' ( B1 ) 



where P„[(0,0) — > (x,t)] is the probability that the line 
(0, 0) — ♦ (x, t) cuts precisely n domain walls, k is the 
number of upsteps, 6(m) is the discrete delta function 
modulo 3, and A/i = fi' — /i. The probability that we cut 
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the path of n particles when moving from (0,0) to (x,t) 
can be computed in terms of cutting none: 



P„((0,0)-> (x,t)) 



M 



P M ~ n (l-P) n , (B2) 



where M is the total number of particles and 



P = P(will not cut v'th) = 



L/2 dx„ 



L/2 



L 



dv v P(v v ) 



dv v P{v v ) 
(B3) 



which can be calculated to get 
P(will not cut v'th) = 1 - ^~ 



—=e " + u crf{u) 



(B4) 



with u the dimcnsionlcss velocity introduced earlier, u 
x/i. 

Using then the representation of the delta function 



5(m) = Re 



1 

3 ' 



all summations above can be evaluated to yield 



(B5) 



(B6) 



leading to the relaxation function given in the main text, 
Eq. 121. 



2. case: Paramagnetic side 

Two differences arise when we use the diagonal <S- 
matrix, Eq. {TTJ. The first one is that the sequence of 
the colors now does change in course of the scattering 
process, while the quantum number of a given line does 
not. The second change is that the number of the (-1) 
factors picked up by the wave function now depends on 
the quantum numbers of two colliding particles. This sec- 
ond change implies that, after averaging over the color of 
the other quasiparticlcs, one gets identically if the par- 
ticle created by the operator P M collides with any other 
particle. Thus the relaxation function is simply the prob- 
ability that the created particle propagating from point 
(0, 0) to (x, t) collides with no other particle. 

The probability that the z/'th particle will not collide 
with the injected one is given by Eq. (|B4(1 . 

Thus the probability that the injected particle does 
not collide with any of the already existing quasiparticles 
simply reads 

^TgrableOMl = Y[ P(will not Cut v'th) = 



= 1- 



tG(u) 
M 



M 



,-tG(u) 



(B7) 



resulting in the relaxation function given by Eq (|23|) . 



APPENDIX C: RELATION BETWEEN THE 
EFFETIVE HAMILTONIAN AND THE 
SCATTERING MATRIX 

Already simple rcnormalization group arguments sug- 
gest that the 5-matrices above describe the scattering of 
quasiparticles with vanishing momenta for any coupling 
g 7^ g c and T <C A: We know very well that the cou- 
pling g is relevant for g > g c and scales to g — ► oo under 
the RG flow , while for g < g c it is irrelevant and scales 
to g — ► 0. The asymptotic (long wavelength) dynamical 
properties, however, remain invariant under the renor- 
malization group, and therefore the scattering matrix of 
quasiparticlcs with vanishing momenta obtained by sim- 
ply performing perturbation theory in 1/g (or g) around 
the trivial fixed point Hamiltonians must coincide with 
the exact scattering in both phases. 

One can refine the rcnormalization group argument 
above to show that, apart from some very special points, 
the iS-matrix of the quantum Potts model on a lattice 
should always take on the simple form, Eq. HOfl for van- 
ishing quasiparticle momenta as follows. Let us suppose 
that we have a high energy cut-off A, larger than the gap, 
A > A. Performing a rcnormalization group transforma- 
tion down to a length scale b ^> A -1 , we obtain a local 
Hamiltonian for the elementary excitations ('dressed' lo- 
cal flips and kinks for g > g c and g < g ci respectively), 
that we can simply construct based upon symmetry con- 
siderations and power counting: 



H. 



off 



E- 



d 



2A dx % 2 
irrelevant terms 



E 



6(xi - Xj) 



(CI) 



where Xi denotes the coordinate of the i'th quasiparticle. 
Here the second interaction term is always relevant, and 
it is this term that determines the structure of the 5- 
matrix for vanishing momenta. For the sake of simplicity, 
let us restrict ourselves to the most interesting Q = 3 
case. Then A is just a chirality label, A = ±, and the 

u , u 2 



interaction matrix u x '\ J . can be characterized through 



three dimensionless parameters, u\ = u +H 



and M3 



"++ 



These parameters 



are dimensionless functions of A b and A/ A and can be 
written as 



Ab, 



A 



Q 



1,2,3 



(C2) 



In the b — > oo (large wavelength) limit these must scale 
to constants that only depend on the ratio A/A, 



lim u 

b — >oo 



(C3) 



These three numbers determine then the asymptotic form 
of the 5-matrix. It is then a simple matter to show that 
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the (S-matrix generated by this effective interaction al- 
ways assumes the form Eq. llUl a^ excepting the special 
case when 

A = u\ . (C4) 

In this latter case the interaction is singular (vanishes 
in one of the scattering channels) and one obtains an 
iS-matrix that is diagonal for vanishing quasiparticle mo- 
menta, as given by Eq. 111|) . This S- matrix coincides 
with the one that one obtains by requiring integrability, 
i.e. by requring that the scattering matrices satisfy the 
Yang-Baxter relations^ Thus Eq. I|C4(1 can be viewed 
as an integrability condiction. However, the Q = 3 state 
Potts model is not integrablc away from the critical point. 

Although there is no obvious reason for that, never- 
theless, it might be possible that Eq. I|C4() is satisfied 
if one removes the cut-off from the theory, i.e. if one 
takes the limite A/A — * oo. Indeed, the scattering ma- 
trix (|llf> is only obtained in a continuum field theory ap- 



proach, when one removes the cut-off and also assumes 
integrabilityi 31 i 44 i 45 . 

However, in a theory with a cut-off there is not a single 
reason why the integrability condition, Eq. I|C4|) , should 
be satisfied. In fact, our pcrturbative expansion just 
proves that Eq. I|C4|I is not satisfied in the large g and 
small g limits. 

These arguments very convincingly support that the 
5-matrix of the quantum Potts model on a lattice al- 
ways takes on the the form (|l(Jf) in the limit of vanishing 
momenta, and the low energy fixed point theory is simply 
generically not integrablem^ It is an interesting question 
if there exists still a cross-over regime for A <C A, i.e., 
if there is some regime where for intermediate energies 
A < c 2 q 2 /A ~ T > T* , the diagonal 5-matrix is ade- 
quate, but the discussion of this requires extensive nu- 
merical studies and is certainly beyond the scope of the 
present paper 4^ 
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